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I. INTRODUCTION 


The underwater echo-sounder is used mainly for sounding purposes and has 
become an essential part of the equipment used by hydrographers. But sounding 1s not 
the only role for echo-sounders; it may also serve as a fish-finder for fishermen, sub- 
bottom profiler for marine geologists, etc. 

To provide nautical charts containing information about the bottom type, 
acoustic techniques can provide a rapid remote method for the examination of the sea- 
bed and sub-bottom. It has been and continues to be the aim of many workers to 
apply acoustic methods to the quantitative classification of sedimentary sea-bottoms. 

A sound beam projected onto to the sea-bottom interacts with the water- 
sediment interface, the nature of its echo will depend on the bottom type and 
roughness. For’plane fluid reflectors, the amplitude of an acoustic wave echo is 
governed by the Rayleigh reflection coefficient (Bezdek, 1973). The coefficient 1s an 
explicit function of the properties of the two mediums, 1.e. their densities, sound speeds, 
and the angle of incidence. 

However, if the boundary is not plane, instead acting as a scattering surface, 
scattering will affect the value of the apparent reflection coefficient; the degree that it 
does depends on the scale of the boundary irregularities relative to the acoustic 
wavelength (Huang, 1971). In this situation the result is that some energy will be 
scattered backward from the interface and some will penetrate into the bottom, where 
It is both absorbed and scattered. The backscattered energy from within the bottom is 
- manifested as. bottom reverberation. This is illustrated in Figure 1.1 which shows 
sound impinging on a rough bottom where reflection occurs at the interface before and 
after it is scattered by volume inhomogeneities in the sediment. 

Since different sediments have different grain sizes, they have different surface 
and volume reverberation characteristics. How to use this simple but important fact to 


classify the materials forming the surface of a sea-bed is the subject of this research. 
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Figure 1.1 Bottom acoustic scattering mechanisms. 


Il. BACKGROUND AND LITERATURE REVIEW 


To provide an understanding of this experiment some background information 


relevant to the interaction of sound with the water column and sea-bed is provided. 


A. REVERBERATION 

In an empty room or in a mountain valley it is not unusual for a sound to persist 
even though its source has been cut off. This reverberation usually decreases rapidly in 
intensity with time. When distinct echoes are not heard, the reverberation may be 
considered as a blending of a large number of echoes. 

Similarly, when a sound beam is projected into the ocean from a sound projector, 
a persistence of sound is observed either when the projector serves also as a 
hydrophone (monostatic case), or when a separate hydrophone is used (bistatic case). 
But since the ocean is really not a.closed volume like a room, and the medium is also 
different, the explanation of ocean reverberation must differ from that for room 
reverberation. 

In the ocean: the roughness of the surface and bottom, air bubbles, suspended 
solid matter, and organic matter such as plankton and fish cause inhomogeneities. 
These “targets” scatter sound back to the hydrophone as a large number of weak 
echoes which are observed as reverberation or backscattering (Eyring, 1948). 

In an ocean environment the reverberation may be one of three types according 
to where the scatterers are located. One type of scattering occurs in the volume of the 
sea, or in the volume of the sea sediment, and produces volume reverberation. Sea- 
surface reverberation is produced by scatterers located on or near the the sea surface, _ 
and bottom reverberation originates at scatterers on or near the sea bottom. The last 
two types of reverberation may be considered together as surface reverberation, since a 


two-dimensional distribution of scatterers is involved (Urick, 1983). 


B. HISTORICAL PERSPECTIVE 

The investigation of sound scattering by the ocean bottom began during the 
Second World War and was carried out both in shallow and in deep water (Bunchuk 
and Zhitkovsku, 1980). The earliest attempts to use backscattered sound to classify 
sediment types were at both normal angles of incidence and at different grazing angles. 
Frequencies used were mostly below 200 kHZ (Gurcan, Creasey and Gazey, 1981). 
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Descriptions of the backscattering fall into two categories: One approach is to 
attribute the backscattering to surface roughness, while the other is to consider that the 
bottom backscattering is due to volume backscattering ( Eckart, 1953; Horton and 
Muir, 1967). 

A different approach was followed by Pace and Ceen (1982) who showed on the 
basis of the time spreading of rectified returns that sea-bottoms of mud, sand, fine 
gravel, coarse gravel may be distinguished one from another. Gurcan, Creasey and 
Gazey (1981) found the power spectrum of backscattered sound to exhibit 
harmonically-related resonance peaks. From their positions and spacings these 
resonances can be used for classification purposes. 

Other important work was done by Breslau (1967), McLeroy (1972), Dodds 
(1984), and Meng and Dinghua (1984). Breslau was concerned with interference based 
on the reflection coefficient of the bottom; McLeroy and Meng and Dinghua were 
concerned with extracting information from the echo waveform. Dodds company, 
Huntec, Ltd., has developed the Deep Tow System which uses a broadband acoustic 
source (1-10 kHZ) to remotely classify the bottom with some success. Recent work 
done by Diaz (1986), Chang (1986), and Facada (1987), used high-frequency (181 kKHZ) 
tone-burst signals at normal incidence to show that the volume reverberation of 
sediment is proportional to grain sizes. 
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It. MODEL DEVELOPMENT 


Since sound waves interact with the water-sediment interface, the backscattering 
energy from the interface is dependent not only on the surface roughness but on the 
sediment volume scattering as well. The theoretical considerations of this experiment 
therefore should include both surface and volume reverberation from the sediment. 
For this reason some basic surface and volume scattering theory will be discussed, and 


finally, a simple model will be developed. 


A. SCATTERING STRENGTH AND REVERBERATION LEVEL 
1. Scattering Strength 
The physical property upon which reverberation depends is called the surface 
or volume scattering strength. It is defined to be 


S.y = 10 log (Iscat/linc) (3.1) 


In this parameter the unit is decibel, I,.., is the intensity of the sound 
scattered by a unit area or volume refer to a unit distance, and I; is the incident 
plane wave intensity. If the direction of scattering is back to the source, the quantity S. 
is termed backscattering strength. The scattering strength for volume and surface 
reverberation are illustrated in Figures 3.la and 3.1b (Urick, 1983). It is customary to 
‘use the backscattering cross section of a unit area or volume as the basic parameter 


involving scattering; it is given the symbols m, and m, and is the ratio of scattered 


Vv 
power, refered to unit distance, to the intensity incident on a unit volume or area. 


Thus, 


I Mg y/ area or volume (2) 


scat 


1s. 


If it is assumed that the scattering is uniformly distributed over a sphere, or a 


hemisphere, Equations (3.1) and (3.2) give 


” 
I 


y = 10 log (m,/47) . (3.3) 


Y 
I 


10 log (m,/27) (3.4) 


2. Reverberation Level And Assumptions 

In the sonar equation, the term RL refers to the “equivalent plane-wave 
reverberation level”. This is the level of the axially incident plane wave which produces 
the same hydrophone voltage across the hydrophone terminals as that produced by the 
received reverberation. 

To find an explicit expression for RL, some simplifying assumptions are 
needed: 

1. Straight-line propagation paths, with all sources of attenuation other than 
spherical spreading neglected. 

2. A random, homogeneous distribution of scatterers throughout the area or 
volume producing reverberation at any one instant of time. 

3. A density of scatterers so large that a large number of scatterers occur in an 
elemental volume dV or dA. | 

4. A pulse length short enough for propagation effects over the range 
extension of the elemental volume or area to be neglected. | 

_5. An absence of multiple scattering, that is, the reverberation produced bv 


reverberation is negligible. 


B. SURFACE SCATTERING THEORY 

Surface scattering means scattering produced by scatterers distributed over a 
nearly plane surface rather than throughout a volume. In terms of decibels the 
reverberation level can be expressed as 


RL, = 10 log [(I,/r*)u, f b(8,9) b’(0,9)dA] (3.5) 
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Figure 3.1 scattering strength for volume and surface scattering. 


17 


Where 
dA = elemental area of the scattering surface. 
I, = axial intensity at unit distance from source. 
r = distance between transducer and scattering surface. 
10 log u, = S, is the scattering strength for surface 
reverberation. 
b(8,9) = beam pattern function of the projector. 


b’(8,9) = beam pattern function of the hydrophone. 
® = vertical beam angle (normally small) 
? 


b(8,@) is the beam pattern function which acts to define the projector’s major 
lobe. On the beam center b(8,@) = 1, but it drops to 0.707 at the 1/2 power point, i.e. 
I, is reduced by 3 dB. In monostatic sonar the same transducer is used both for 


horizontal beam angle. 


projection and reception, and we can assume that transmitted and received patterns are 
the same; but in practice they need not be. 

If dA is taken as a portion of a circular annulus in the plane of the scatterers, 
according to Figure. 3.2 we may then write 


dA = (ct/2)rdO os (3.6) 
Then reverberation level becomes 
RL, = 10 log (Ig/r*u,(ct/2)r {o(8,¢)b’(8,9)d© : Ba 


For near grazing incidence, the axis of the beam is only slightly inclined toward 
the scattering surface, so that 


J b(0,9)b'(0,9)do = f b2(0,e)do = ® (3.8) 
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Figure 3.2 Elemental area of the scattering surface. 
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Where @ is the effective horizontal beam width of the transmitting/receiving 
transducer, therefore 


RL, = 10 log [(1,/rPu,(ct/2)r®] (3.9) 


= SL-2TL tes, + 10dcees 


Where 
10 log I, = SL 
10 log 1/r4 = - 2TL 
10 log u, = S, 


(ct/2)r® = A 


C.. VOLUME SCATTERING THEORY 

Due to inhomogeneities within the ocean volume, sound is scattered in all 
directions, and the portion of scattering which returns to transducer is known as — 
backscattering or reverberation. Sound which has penetrated into sediment and is 


scattered back is also considered as volume reverberation. Volume scattering is 
ulustrated in Figure 3.2a and 3.2b. | 


The volume reverberation can be expressed as 


RL, = 10 log [(I,/r*)uy J 6(8,9)b’(6,)dV)] (3.10) 


Individual terms have the same definition as used for surface reverberation except 
for the exchange of S,, for S, and dV for dS. Here 10 log u, is the backscattering 
strength for volume reverberation and is denoted by the symbol S,. The elemental 
volumes according to Figure 3.3 are 


dV = r(ct/2)dQ | (3.11) 
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Figure 3.3 Geometric view of volume scattering. 
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Where T is pulse length, c is the velocity of sound, and dQ 1s the elemental solid angle 


subtended by dV at the source. The volume reverberation 1s 


RL, = 10 log [(I,/r4) r?(ct/2)S, J bb’dQ} (3.12) 


Let the actual two-way beam pattern bb’ be replaced by an ideal pattern as 
shown in Figure 3.4. Then the width of this ideal beam is given by 


fob'dQ = J 1x1d2@=¥ (alg) 


In terms of ‘P the reverberation level becomes 


RL, = 10 log [(Ip/r*)uy(ct/2)¥r7] (3.14) 


= SL-2TL +S, + 10logV 


where 
10 log (Giese 
10 log I/r4 = - 2TL 
lO leg qu S 


(ct/2)¥r2 = V 


In summary the reverberation levels for surface and volume are - 


RL, = SL-40logr +$,+10logV 9 V = (ct/2)¥r? 
RL, = SL-40logr +S, + l0logA A = (ct/2)@r 
The most complicated terms in these two equations are the scattering strength S, and 


Sy in which 


S77 Legs S, = 10 log u, 
The quantities of u, and u, can depend on the grazing angle of the incident wave 
relative to the bottom, the frequency, the scattering angle, and the surface relief and 
structure of the bottom (Buhchuk and Zhitkovsku, 1980). 
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Figure 3.4 The elemental volume for volume scattering. 
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Figure 3.5 Actual and equivalent beam patterns. 
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D. MODEL OF ACOUSTIC BACKSCATTERING 


The model of Sanders and Coppens (1986) was used for this experiment. It is as 
follows 


The square pulse used to excite the transmitting transducer can be described by 
the unit step function, H. 


V(t) = V,[H(t) - H(t-t)] (3.15) 


where V(t) is the voltage at time t, V, is the maximum voltage, and t is the pulse 
length. 


0 t<0 
mo) = 

1 t20 

0 t-t< 0 
H(t-t) = | 

1 t-7 20 


The square electrical pulse is transformed by the transducer into an acoustic 
pulse with an exponential rise, an approximately flattened top, and an exponential 
decay. This shape of waveform can be described as 


f(t) = (1-e"@5 H(t) ; (1-e°(t 7) H(-2) | (3.16) 


where t is the time measured from the beginning of the pulse and @ is the rise and 
decay constant (assumed to be the same). 


For a transducer located in the water and sending out a sound pulse to a smooth 


sediment surface, the reflected pulse received by same transducer will have an intensity 


2 


1 = (f T.2 5, PS 1 (t-26/cydb(r+&)* + Rikt)/r? yr? (3.17) 


where 


qT; 


intensity transmission coefficient 


s, = volume scattering coefficient 

8 = absorption coefficient in sediment 
€ = affected sediment depth 

R; = intensity reflection coefficient 


r = distance between transducer and sediment 


From (eqn 3.17), the intensity contributed by scattering is 

1 = § 7.2 se PS 1, mt-26/c)db (r+ &)2/r2 (3.18) 
and the intensity contributed by the reflection is 

L. = RI, ftt)/r4 (3.19) 


_ Ifwe assumer > > 6, then 


1. = (Ts, Ly/r9) § ec PS mt-28/c) db (3.20) 
where 
f(t-2E/c) = (1-e78(t25/¢))1(t-26 /c)-( 1-076 7-28/ )yHH(t-2-2E/c) (3.21) 
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The lower limit of & is zero; the upper limit is ct/2 for H(t-26/c) and c(t-t)/2 for H(t- 
t-26/c). Here only the decay portion, t27, of interest. After integration we have 


1 = A {(1-e7Bt/2)/p . (e-CBU/2_e-@tyy(-B + 2a/c) (3.22) 


+ (e°B(ET/2 1/8 + (e°CBLt-8)/2_¢-A(t-T)) / (8 + 20/c)] 


where 


A= T.* sy 1, /e° (3.23) 
with f(t) from Equation 3.16, I. becomes 

I. = 1, R; (e%* -1)e@tr4 (3.24) 
Adding - components I, and I, gives the total intensity | by the transducer: 


I = N{T,2s,/Ple"Pt/2(eCB2/2_1)) + 7.25, /(-B + 2a/c) (3.25) 


fe" CBt/2/.cBt/2_ 1)+ e At |-e7)] a Re. De 


where N = | ait" 


This formula looks quite complex, but all the variables in this expression are 


Known except s, . For this reason, trying different s, values to get one to fit the 


V 
experimental data for different kinds of sediments is an important goal of this research. 
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IV. EXPERIMENTAL DESIGN AND PROCEDURES 


A. EQUIPMENT AND ACCESSORIES OF EXPERIMENT 
1. Material Selection 
In this experiment, five different sediments were used: Monterey fine sand 
#30, Monterey Aquarium #2, aggregate, and two sizes of glass beads. The sizes of the 
materials varies from 0.2 mm to 5.3 mm (Facada, 1987) and are listed in Table 1. 
Since the natural sediments used are inhomogeneous in composition and shape, a 
“sediment” composed of homogeneous particles (spherical glass beads) was chosen for 
some experiments. : | 
The various materials were added to water and left in the test tank for two or 
three days to allow the water to completely saturate the sediment and to let air bubbles 
dissolve. To help remove air bubbles, a water jet or hand stirring was used until no 
further bubbles could be seen rising from the sediment. Bleach was added to the water 
to control biologic growth. 
2. Experimental Tank . , 
Two glass tanks, each measuring 70 cm X 70 cm X 60cm, were filled with 24 
cm of fine sand, aquarium sand, or aggregate sand,~and 36 cm of water was added 
above the sand. Because of the limited volume of glass beads available, a plastic 
bucket was placed in the tank to contain them. 
3. Electronic Equipment | 
A block diagram showing the equipment configuration is shown in Figure 4.1. 
The output from a General Radio Model 1310 oscillator with a frequency range of a 2 
Hz to 2 MHz was simultaneously fed into a Hewlett- Packard HP 5233L frequency 
counter and a General Radio Type 396-A tone burst generator which was used to 
generate either 8- or 16-cycle pulses. The signal output passed through a Hewlett- 
Packard HP 467-A power amplifier, into a Datasonic Transmit/Receive (T/R) switch, 
then to the transducer. The output signal was also used to tngger a Tektronix Type 
RM 503 analog oscilloscope and a Nicolet Model 3091 digital oscilloscope. 
The signal return from the sediment was received by the transducer, amplified 
20 dB by a Hewlett-Packard 465-A pre-amplifier, and then passed through a Spencer 
Kennedy, Inc., Model 302 variable electronic filter which serves as a high-pass filter to 
eliminate low frequency noise before going to the oscilloscope. 


28 


TABLE 1 
oe wile ROPERTIES OF SEDIMENT 


FINE AGGREGATE AQUARIUM COARSE FINE 
SAND SAND SAND GLASS BEADS GLASS BEADS 
a Ges 5.3 1972 or | 0.16 
( mm) 
p 2.69 2.64 2.59 2.50 2248 
(g/em? ) 
Pee, 1-98 1 4917 22 108 iL Sle 
(g/cm? ) 
e 1680 = 1555 1588 1822 1878 
(m/s) 
cae 0.36 0.36 0.36 0.53 0.48 
i, 0.198 0.52 0.38 O41 0. 42 


Three different transducers were_used in the experiment. A Naval Research 
Laboratory F-41 circular-piston type had a 8.8 cm diameter active face and a thickness 
of 4.4 cm. Its resonance frequency was 185 kHz (Diaz, 1986). The second, a Ross 
transducer, was borrowed from the Acoustic Division of the Atlantic Meteorological 
and Oceanographic Laboratory (AMOL), Miami, of the National Oceanic and 
Atmospheric Administration (NOAA). This seven-element transducer has a frequency 
range from 20 kHz to 150 kHz with resonance frequencies at 62 kHz and 108 kHz; the 
-3 dB beamwidth at 62 kHz is 10 degrees, and at 108 kHz it is 5.5 degrees 
(AMOL/Ocean Acoustics Division, 1979). The third transducer, a Celesco Industries 
Type LC10 hydrophone, was used to measure the sound speed in the sediments. The 
LC10 is a small (0.97-cm diameter by 2.87-cm long) cylinder with a receiving range of 
0.1 to 120 kHz. 
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Figure 4.1 Schematic of the equipment configuration. 
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The waveform of the return signal was displayed on the digital oscilloscope 
and sent as a stream of characters to a HP-86 computer. To reduce manual operations 
a HP-3421A Data Acquisition and Control Unit was used to receive commands from 


the computer and then automatically control the digital oscilloscope. 


B. MATERIAL PROPERTY MEASUREMENT 
The physical properties of the sediments were measured for later use in 
computations. 
1. Density 
The water-saturated and dry densities of the different sediments were measured 
separately. The dry density was measured by filling a 100-ml graduated cylinder and 
weighting. After subtracting the weight of the empty cylinder to find the net weight of 
material, division by the volume gives the density. Water-saturated density was 
measured using the same method. The density measurement results are summarized in 
Table 4.1. 
2. Speed of sound 
The speed of sound was measured by using three LC-10 transducers inside the 
tank, one acting as source, the other two as receivers. The three LC-10s were inserted 
at the same -depth in the sediments. One LC-10 was placed 15 cm away from the 
source, and the second LC-10 was moved along the straight line defined by the source 
and first receiver. The two signals shown on the digital oscilloscope (Figure. 4.2) allow 
measurement of the travel time over the known distance between two LC-10 receivers. 
From this the speed of sound in the sediment can be calculated. 
3. Reflection coefficient | 
The normal-incidence pressure reflection and transmission coefficients were 
measured to check the computed value. The ratio of maximum measured voltage for 
signals returned from the sediment and from the water/air interface at same distance 


gives the reflection coefficient, 
eee Vie) Va - (4.1) 


where 
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= Ae 


Figure 4.2 Two echos received on the digital oscilloscope. 
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Vy = maximum voltage of signal return from water/air interface. 
V, = maximum voltage of signal return from sediment. 


Since the densities and sound speeds are known, the normal reflection 


coefficients of the different sediments can also be computed from 
fe (89 Ey) (to + Ty) ) (4.2) 


where 


rj = Py Cj is the specific acoustic impedance of water. 
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PC is the specific acoustic impendance of sediment. 


The measured and computed reflection coefficients are listed in Table 1. 


C. EXPERIMENTAL DESIGN CONSIDERATIONS 
1. Fundamental concept of experiment 

The original hypothesis of this experiment was that the different sediments 
must have different volume reverberation characteristics and that there should exist a 
relationship between the sediment properties and their volume reverberations. 

Let a transducer located in the water send out a sound pulse to the water/air 
interface. -If the interface is smooth, the received echo waveform is purely reflection. 
If the interface isnot smooth, then the amplitude of the received waveform is reduced 
because the sound energy scattered at the rough surface. Suppose the scattering is first 
order and the surface has the same roughness everywhere; then the echo is a 
combination of reflected and scattered energy, and its amplitude is reduced in relation 
to the purely reflected echo amplitude. In Figure 4.3(a) the solid line represents a 
purely reflected waveform, and the dashed line represents the combined reflected and 
scattered waveform. 

If the echo is returned from a water/sediment interface instead of the water/air 


interface, the echo amplitude in the tail should be increased because additional energy 


be 


is scattered back from within the volume of sediment. In Figure 4.3(b) the dotted line 
represents the echo from the sediment, while the solid and dashed lines represent the 
purely-reflected water/air interface amplitude and the combined reflected and scattered 
amplitudes from the water/air interface. 
2. Experimental Considerations 
Some necessary considerations are required to ensure the success of this 
experiment. Some of these effects were studied by Diaz (1986) and Chang (1986). 
a. Effects of the transmitting electronics 
The the electrical signal from the tone burst generator applied to the 
transducer showed that the signal envelope was a clean square wave. This 
demonstrated that the transmitting electronics had a negligible effect on the received 
signal. 
b. Effects Of Sampling Rate 
The fastest sampling rate of the digital oscilloscope was 1 ps. and is not 
enough to give a precise reproduction of the echo waveform at the high frequencies 
used for this experiment. However, observation showed that there was natural “jitter” 
in the triggering of the scope, so that for consecutive pulses, the scope would sample 
different parts of the waveform. To obtain an accurate waveform, a large number of 
waveforms were sampled and averaged. If sufficient “jitter” were present, the voltage 
measured in a given sampling bin should vary between the maximum and minimum 
values in the vicinity of the corresponding time. Since the sine wave varied between 
negative and positive, these voltage values were squared before averaging in order to 
get an accurate shape of the waveform envelope. | 
To determine the minimum number of samples required for good 
. reproduction of the shape of the waveform envelope, the water/air interface was 
measured with different number of samples by Diaz (1986) and Chang (1986). Their 
results showed that 50 samples for each data set were enough to give a true echo 
waveform. 
c. Effects Of Water] Air Interface Interference 
In an ideal situation the transducer should be far away from all unwanted 
surfaces to ayoid interference. If the transducer is located close to the water surface 
and projects a pulse to the sediment surface below it, most of the energy will be sent to 
the sediment surface, but nevertheless some energy will be transmitted to the water 
surface and then reflected back toward the sediment. Restricted by the size of the test 
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(b) ECHO WAVEFORM FOR SEDIMENT SURFACE 


Figure 4.3 Theoretical echo waveform from water and sediment. 
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tank, the transducer had to be located close to the water surface. To find out the 
magnitude of this interference, the water surface was agitated by hand to alter the 
amount of reflection from the water surface. The observed waveform was not changed 
during the ruffling. It was concluded that this interference is negligible. 

d. Effects Of Bottom Reflection Interference 

To find out if the sediment in the test tank was thick enough to avoid 

reflection from the bottom of the tank, a simple computation is sufficient. For a 
16-cycle signal with a frequency of 181 kHz, the pulse duration is 88 ws, and the decay 
time approximately 100 jis; the echo from the smooth water/air interface has the same 
decay time. For sound pulses in the sediment the decay time should be longer because 
of volume reverberation. For the fine sand the decay time is 152 ts, for a total time 
interval of 240 tts. With a sound speed ‘of 1680 m/s for the fine sand, and a sediment 
depth of 25 cm, the echo from the bottom of the tank comes back to the surface of 
sediment 286 pts after the incident pulse begins reflecting from the surface of the 
sediment, too late to interfere with the reflected pulse. For the other sediments the 
depth was deep enough to prevent any bottom reflection interference. 

e. Effects Of Side Lobe 

- Since the test tank was too small to avoid the interference, if any, caused 
by the side lobes, a test was necessary to find out if this is important. The half- 
beamwidth (-3 dB) for the F-41 transducer was measured to be 10.0° (Chang, 1986); 
for the Ross transducer the half beamwidth at 108 kHz it was 5.5° and at 62 kHz it 
was 10.0° (AMOL, 1979). For a 30-cm distance from transducer, the sediment surface 
insonified by the major lobe is a circular area of radius about 5.3 cm. Since the test 
tank dimensions were 70 cm X 70 cm X 60 cm, if the transducer is located at the 
center of tank; the major lobe will not strike the sides of the tank. To‘test this, Diaz 
(1986) and Chang (1986) disturbed the surface of the sediment outside the center 30 cm 
x 30 cm square area. Waveforms returning from the sediment were compared before 
and after the disturbance. They were indentical, indicating that the effect of the side 
lobe was negligible. For the Ross transducer the same measurement was: made, 
resulting in the same conclusion. | 

i. Effects Of Sediment Inhomogeneity 

The’ experiments of Diaz (1986) and Chang (1986) were bothered by a 

drastic change in the echo waveforms whenever the transducer was moved parallel to 


the smoothed aggregate or fine sand sediment surfaces. It was assumed that the 
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waveform fluctuations were caused by inhomogeneities in the sediment. Facada (1987) 
tried to use more homogeneous glass beads as sediment to prove this assumption, but 
the results did not satisfy him; the inconsistencies between echo waveforms from 
different pulses still existed. He concluded that the variations in the shape of the pulse 
tails increase as the grain size of the sediment increases and that changes in the echo 
waveforms resulted not from volume reverberation but from surface scattering. 

To avoid particulate scattering the wave lengths should be at least three 
times longer than the grain size (Bradshw, 1981). To meet this criterion three different 
frequencies were used. For fine sand, since the largest grain size is 0.6 mm, a 
wavelength equal to 9.3 mm (181 kHz) readily qualifies. Two other frequencies, 108 
kHz and 62 kHz, are applied to aquarium and aggregate sand, respectively. 

g. Effects Of Gas Bubbles In The Sediment 

Since we are only interested in the decay constant of the sediment itself, gas 
bubbles trapped inside the sediment should be removed before the experiment. To help 
remove gas bubbles a water jet was used to stir the sediments before each set of data 
measurements. | 

h. Effects of Gas Bubbles on the Active Face of the Transducer 

When a cloud of gas bubbles collects on the active face of the transducer, 
the sound waves will be affected and the amplitude of the echo will be reduced. This 
may be seen on the oscilloscope. It was very important to remove these bubbles from 


time to time during the experiment by means of a small brush. 


D. PROGRAMMING | 
The computer programs used in this experiment were prepared by Chang (1986). 
To acquire and analyze the data three programs are available. 
1. Program ’THESIS3” | 
This program converts the waveform data from the Nicolet 3091 Digital 
Oscilloscope into voltage and time value. The Nicolet is controled by the HP-3041 
Data Acquisition and Control Unit. The program squares the voltage and then takes 
the average of 50 sets of readings. To insure the data are not lost, the data are stored 
on disk after each 10 sets are acquired. | 
2. Program “MAXIMUM” 
This program reads the squared voltage data and the corresponding time data 
stored on disk and then determines local maxima. The square-roots of the maxima 


were taken and stored on the disk in a separate data file. 
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3. Program ”“PLOT3” 
This program plots the envelope of the echo waveforms. The voltage is 
converted to log values and plotted vs. time. A least-square method 1s used to find the 
best fitting straight line, which is then plotted. The slope of this line and its correlation 


coefficient are then printed. 
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V. RESULTS AND DISCUSSIONS 


One of the significant features of this experiment is the measurement of volume 
reverberation for different sediments. The theoretical echo waveforms for water and 
sediment can be seen in Figures 4.4(a) and (b). Actual echo waveforms from water and 
sediment are shown for comparison in Figure 5,1. 

Since the observed volume reverberation in sediments is quite sensitive to 
inhomogeneities in both bottom surface and composition, the echo waveforms vary not 
only from sediment to sediment but also between individual measurements on the same 
sediment. For this reason, 50 samples for each data set and 10 data sets for each 
sediment type were taken to get the average shapes of the echo waveforms. 

The slopes of the trailing edges of the echos for the various sediments were 
measured, and a relationship between grain size and the echo decay slopes at different 
frequencies was determined. Finally, a simple model of backscattering was used to 


determine the scattering coefficient for each sediment. 


A. DECAY OF VOLUME REVERBERATION 
1. Fine Sand | 

Ten sets measurements were made at 181 kHz for fine sand (ka = 0.11) and 
the individual echo decay slopes found along with their correlation coefficients (Table 
a 

The decay slopes for fine sand varied from -0.0310 to -0.0378 Np/pts. The 
average slope of-all ten sets was -0.0344 Np/ps. Compared to the decay constant for 
purely-reflected (water/air interface) pulses, a = -0.0377 Np/ps, volume reverberation 
- 1s seen to be present. 
Since the fine sand was the most homogeneous sediment used and its surface 


was smooth, the variation in its decay constant was smaller than for the others. 
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Figure 5.1 The volume reverberation of sediment. 
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TADBEE 2 
FINE SAND INDIVIDUAL SLOPES AT 181 KHZ 


DATA SLOPE x 1074 CORRELATION 
SET (Np/us) COEFFICIENT 
DATA10 =e), 06 -0. 9648 
DATA11 -3.38 -0. 9548 
DATA12 =3 5110 -0. 9605 
DATA13 eae -0. 9817 
DATA14 = oemiel -0. 9814 
DATA15 ~3, 62 -0. 9854 

- DATA16 ~3. 37 -0.9654 
DATA17 ~3.57 , -0. 9545 
DATA18 ~3, 44 -0. 9823 
DATA19 =30 9 =0. 99726 


MEAN SLOPE = -3.44 x 1074 Np/ps 
¢ = 0.22 x 1074 Np/ps 


2. Aquarium Sand 

Two different frequencies, 62 kHz and 108 KHz, were used for the aquarium 
sand experiments. Individual and averaged slopes are listed in Tables 3 and 4. 

At 62 kHz an 8-cycle of pulse is 129 pts long, and its ka value is about 0.21. 
At 108 KHz an 8-cycle pulse has a duration of 74 ts and a ka value of about 0.37. 
Both may be considered within the Rayleigh regime. 

The decay slope of this sediment varies from -0.0035 Np/Ms to -0.0085 Np/us 
at 62 kHz and have average slope of -0.0059 Np/ps. At 108 kHz the slope varies from 
-0.0152 Np/ps to -0.0208 Np/ps with the average being -0.0184 Np/ws. 

As expected, the echo waveforms for aquarium sand are less consistent than 
for fine sand. This is probably due to the fact that the: grain sizes are less 


homogeneous and the sediment surface is less uniform. 
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TABLE 3 
AQUARIUM SAND INDIVIDUAL SLOPES AT 62 KHZ 


DATA SLOPE x 1074 CORRELATION 
SET (Np/s) COEFFICIENT 
DATA20 =o cr -0. 6961 
DATA21 -0. 83 -0. 8041 
DATA22 -0.53 -0. 7440 
DATA23 -0. 60 -0. 7720 
DATA2 4 -0. 41 -0. 6950 
DATA25 -0. 66 -0. 7292 
DATA26 =0no5 -0. 6817 
DATA27 -0.85 -0. 6671 
DATA28 -0.52 -0. 6885 
DATA29 -0.50 -0. 733¢ 


MEAN SLOPE = -0.59 X 1074 Np/us 
¢ = 0.16 X 107% Np/ps 


3. Aggregate Sand 

As with aquarium sand, frequencies 62 kHz and 108 kHz were used for the 
aggregate sand experiment. The individual and average slopes are listed in Tables 5 
and 6. . 

At 62 kHz the decay slopes have a maximum of -0.0070 Np/pis and a 
minimum of -0.0034 Np/pts; the average is -0.0051 Np/ps. At 108 kHz the maximum 1s 
-0.0099 Np/pis and the minimum is -0.0032 Np/pts, and the average value is -0.0069 
Np/Hs. 
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TABLE 4 
AQUARIUM SAND INDIVIDUAL SLOPES AT 108 KHZ 


DATA SLOPE x 1074 CORRELATION 
SET (Np/ps) COEFFICIENT 
DATA30 eee? -0. 8949 
DATA31 = eed -0. 8921 
DATA32 GAS : -0. 7076 
DATA33 a7 0 mome77 4 
DATA34 = G -0. 8023 
DATA35 ~1.58 -0. 6813 
DATA36 SoG -0. 8148 
DATA37 5. (ay) -0. 9072 
DATA38 | : -1.73 -0. 8396 
DATA39 =99 -0. 8741 


MEAN SLOPE = -1.84 x 107% Np/us 
¢ = 0.20 <x 107% Np/its 


4. Glass Beads 
Two types of glass beads were chosen for this experiment. Their homogeneity 
in size and shape are the reason why they were used. 

~The average grain sizes of the beads were 1.7 mm and 0.16 mm. Because only 
limited amounts of these glass beads were available, they were put into a plastic bucket 
inside the experimental tank. 181 kHz was used with a pulse length about 88 ws for a 
16-cycle pulse. The ka value for coarse glass beads was about 0.60 and for fine glass 
beads about 0.06. Even though the glass beads were homogeneous in size and shape, 
the individual slopes still varied considerably. Thus, the standard deviation for a total 
of 10 sets of measurements was 0.0023 Np/ps for coarse glass beads and 0.0021 Np/ps 
for fine glass beads. Compared with the value for fine sand of 0.0022 Np/pts there is 
little difference. , 

The echo decay slope for coarse glass beads has a maximum of -0.0358 Np/ps. 
and a minimum of -0.0285 Np/pis, and its average is -0.0338 Np/pts. For fine glass 
beads the maximum is -0.0415 Np/ps, the minimum is -0.0342 Np/ps, and the average 
is -0.0363 Np/ps. All these individual slopes are listed in Tables 7 and 8 
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TABLE 5 
AGGREGATE SAND INDIVIDUAL SLOPES AT 62 KHZ 


DATA SLOPE x 1074 CORRELATION 
SET (Np/ps) COEFFICIENT 
DATA40 -0. 64 -0. 8739 
DATA41 10), Sil -0. 7876 
DATA42 -0. 56 | -0. 7409 
DATA43 -0.70 -0. 6146 
DATA44 -0.54 -0. 7687 
DATA45 =On ce ~0, 2558 
DATA46 Zones -0, 74a 
DATA47 -0. 34 -0. 3158. 


MEAN SLOPE = -0.51 X 107% Np/ps 
G = 0.11 x 1074 Np/us 


TABLE 6 
AGGREGATE SAND INDIVIDUAL SLOPES AT 108 KHZ 


DATA SLOPE x 1074 CORRELATION 
SET _ (Np/ps) COEFFICIENT 
DATASO =omge | -0. 3014 
-DATAS1 =Ou0 -0. 6927 
DATA52 -0. 73 i. -0. 7628 
DATAS3 -0.99 -0. 5083 
DATAS4 -0. 45 -0.1711 
DATASS - -0. 89 ~0. 4688 
DATAS6 -0. 65 -0. 3098 
DATAS7 . ~0.74 -0, 6298 


MEAN SLOPE = -0.69 x 107% Np/ps 
¢ = 0.21 x 107% Np/us 
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B. RELATIONSHIP OF DECAYING SLOPE TO GRAIN SIZE 


NEN / 
COARSE GLASS BEADS INDIVIDUAL SLOPES AT 181 KHZ 


DATA SLOPE x 1074 CORRELATION 
SET (Np/jts) COEFFICIENT 
DATA60 ~3.54 ~0. 8174 
DATA61 Sam 6 -0.9159 
DATA62 -3.56 -0. 9800 
DATA63 aano =Qu oan 
DATA64 ; B2cs ~0.9639 
DATA65 ~ Som ~0.9705 
DATA66 ~3.38 -0. 9610 
DATA67 -3.53 | ~0.9451 
DATA68 oor a ~0. 9800 
DATA69 Some -0. 9869 


MEAN SLOPE = -3.38 x 107% Np/ps. 
G = 0.23 x 107% Np/pts 


Materials of different grain sizes, (fine sand, aquarium sand, aggregate sand, fine 
glass beads and coarse glass beads) were used in this experiment. The grain size and 
average decay slopes are shown in the individual tables. From these results an 
important relationship between volunie- reverberation and sediment grain size can be 
drawn; that is, for a given frequency the amount of volume reverberation increases as 
the sediment grain size increases (Figures 5, 6, and 7). This was expected from the 
suggestion of Clark, Proni, Seem and Tsai (1984). 


45 


TABDBE*s 
FINE GLASS BEADS INDIVIDUAL SLOPES AT 181 KHZ 


DATA SLOPE x 1074 CORRELATION 
SET (Np/ps) COEFFICIENT 
DATA70 ~3. 47 ~0. Sam 
DATA71 -3.62 20.9554 
DATA72 #3, 83 -0. 9042 
DATA73 =3.78 -0. 9209 
DATA74 mes -0. 9904 
DATA75 =3Ne2 -0. 9663 
DATA76 ~353 -0, 99am 
DATA77 -3.66 -0. 9768 
DATA78 =3m2 -0. 9658 
DATA79 | Sona -0. 9580 


MEAN SLOPE = -3.63 X 107% Np/us 
c = 0.21 <x 1074 Np/ps 
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Figure 5.2 Average decaying slope vs. sediment grain size at 62 kHz. 
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Figure 5.3 Average decaying slope vs. sediment grain size at 108 kHz. 
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Figure 5.4 Average decaying slope vs. sediment grain size at 181 kHz. 
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C. RESULTS OF SCATTERING MODEL 
In the model used for this experiment, the total intensity received by the 


transducer 1s 


I = N{T.2s, /Ple"Pt/2(eBt/2_1)] +. 7.25, /(-B + 2a/c) (5.1) 


fe-CBt/2/e¢Bt/2. 1) 4 en At] e2T a R(e®*-1je"Mty 


where N = I iro 


All variables are known in this formula except s,. If the part included inside 


braces 1s called M, then 


[ = (19 /,4)M 


The different sediments have different M values; by assuming a series of 
Scattering coefficients and computing Log(M)!/ 2 for each and normalizing it, the 
model envelopes may be determined. A simple basic program was used, and different 
s,, value are used to find the one for which the slope of the model envelope is equal to 
the average slope of experimental data. (Figure 5.5, 5.6, 5.7 and 5.8). Figure 5.9 
showed the relationship between the scattering coefficient and the ka value. 

In this model, at the beginning of decay the reflection term is dominant and is 
independent of Sy while the scattering term is negligible. This means the model value 
‘will not change in the beginning of decay as the scattering coefficient changes. At the 
tail end of the decay the scattering term is dominant and the decay tails change when 


different scattering coefficients are used. 
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Figure 5.5 Decay slopes of the experiment vs. model for the fine sand. 
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Figure 5.6 Decay slope of experiment vs. model for the aggregate sand. 
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Figure 5.7 Decay slopes of experiment vs. model for the aquarium sand. 
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Figure 5.8 Decay slopes of experiment vs. model for the glass beads. 
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Figure 5.9 Sediments scattering coefficient vs. ka value. 
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I. CONCLUSIONS AND RECOMMENDATION 


The following conclusions can be drawn: 


For a fixed frequency, with increase in sediment grain size there is an increase in 
the volume reverberation of the sediment. 


For a fixed grain size, if frequency is increased, the volume reverberation of the 
sediment decreases. 


Scattering coefficients are roughly proportional to ka; that is, if ka increases, 
the scattering coefficient of the sediment is larger. 


In theory the sediment volume reverberation can be used to classify sediment 
types using of the model to find the best-fitting value of scattering coefficient 
for different kinds of sediment. 


In the real world the ocean environment is so complicated and unpredictable 
that the method suggested may not be practicable. 
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APPENDIX A 
DECAY ENVELOPE GRAPHS FOR FINE SAND SEDIMENT AT 181 KHZ 


The decay envelopes for fine sand sediment at 181 kHz (DATA 10 to DATA 19) 
were showed in this appendix. | 
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Figure A. DATA 10. 
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Figure A.2. DATA II. 
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Figure A.3 DATA 12. 
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Figure A.4 DATA 13. 
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Figure A.6 DATA 1S. 
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Figure A.7 DATA 16. 
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Figure A.8 DATA 17. 
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SLOPE = ~.0344 


TIME (MICRO-SEC) 


Figure A.9 DATA 18. 
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‘Figure A.10 DATA 19. 
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APPENDIX B 


DECAY ENVELOPE GRAPHS FOR AQUARIUM SAND SEDIMENT AT 
62 KHZ 


The decay envelopes for aquarium send sediment’ at 62 kHz (DATA 20 to DATA 
29) were showed in this appendix. 
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Figure B.1 DATA 20. 
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Figure B.2. DATA 21. 
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Figure B.3 DATA 22. 
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Figure B.4 DATA 23. 
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Figure B.5 DATA 24. 
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Figure B.6 DATA 25. 
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Figure B.7 DATA 26. 
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Figure B.£8 DATA 27. 
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Figure B.9 DATA 28. 
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CORRELATION COEFF= -. 7338 
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Figure B.10 DATA 29. 
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APPENDIX C 


DECAY ENVELOPE GRAPHS FOR AQUARIUM SAND SEDIMENT AT 


108 KHZ 


The decay envelopes for aquarium sand sedimént at 108 kHz (DATA 30 to 
DATA 39) were showed in this appendix. 
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Figure C.1_ DATA 30. 
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Figure C.2. DATA 31. 
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Figure C.3. DATA 32. 
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Figure C.4. DATA 33. 
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Figure C.5 DATA 34. 
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Figure C.6 DATA 35. 
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Figure C.7 DATA 36. 
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Figure C.8 DATA 37. 
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Figure C.9 DATA 38. 
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APPENDIX D 
DECAY ENVELOPE GRAPHS FOR AGGREGATE SAND AT 62 KHZ 


The decay envelopes for aggregate sand sediment at 62 kHz (DATA 40 to DATA 
47) were showed in this appendix. : 
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' Figure D.1 DATA 40. 
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Figure D.2. DATA 41. 
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Figure D.3. DATA 42. 
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Figure D.4. DATA 43. 
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Figure D.5 DATA 44. 
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Figure D.6 DATA 45. 
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Figure D.7 DATA 46. 
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Figure D.8. DATA 47. 
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APPENDIX E 
DECAY ENVELOPE FOR AGGREGATE SAND SEDIMENT AT 108 KHZ 


The decay envelopes for aggregate sand sediment at 108 kHz (DATA 50 to 
DATA 57) were showed in this appendix. 
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Figure E.l DATA S50. 
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Figure E.2. DATA 51. 
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Figure E.3 DATA 52. 
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Figure E.4 DATA 33. 
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Figure E.5 DATA 54. 


81 


CORRELATION COEFF= -,. 4688 
SLOPE = -. 0089 


LN CC<VOL TAGE “2>) 71/2) 


2 


TIME (MICRO-SEC) 





Figure E.6 DATA 55. 
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Figure E.7 DATA 56. 
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Figure E.8 DATA 57. 
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APPENDIX F 


DECAY ENVELOPE GRAPHS FOR COARSE GLASS BEADS SEDIMENT 
AT 181 KHZ 


The decay envelopes for coarse glass beads (DATA 60 to DATA 69) were 


showed in this appendix. 
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Figure F.1. DATA 60. 
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Figure F.2. DATA 61. 
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Figure F.3 DATA 62. 
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Figure F.4 DATA 63. 
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Figure F.5 DATA 64. 
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Figure F.6 DATA 65. 
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Figure F.7 DATA 66. 
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Figure F.9 DATA 68. 
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Figure F.10 DATA 69. 
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APPENDIX G 
DECAY ENVELOPE GRAPHS FOR FINE GLASS BEADS AT 181 KHZ 


The decay envelopes for fine glass beads sediment (DATA 70 to DATA 79) were 
showed in this appendix. | 
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Figure G.l DATA 70. 
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Figure G.2. DATA 71. 
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Figure G.3. DATA 72. 


91 


= 
4 

om 
& 
Lu} 
5 Be 
< 
j- 
| 
Oo 
> 
Vv 
ed 
Ww 
Bj 


LN ¢¢(<VOLTAGE “2>) “1/2) 


CORRELATION COEFF= 


SLOPE = -. 0378 


CORRELATION COEFF= 


SLOPE = -.0415 


TIME (MICRO-SEC) 


FigureG.5 DATA 74. 
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Figure G.6 DATA 75. 
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Figure G.7 DATA 76. 
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Figure G.8 DATA 77. 
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Figure G.9 DATA 78. 
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Figure G.10 DATA 79. 
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